Optimizing the fast Rydberg quantum gate 
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The fast phase gate scheme, in which the qubits are atoms confined in sites of an optical lattice, 
and gate operations are mediated by excitation of Rydberg states, was proposed by Jaksch et al. 
Phys. Rev. Lett. 85, 2208 (2000). A potential source of decoherence in this system derives from 
motional heating, which occurs if the ground and Rydberg states of the atom move in different optical 
lattice potentials. We propose to minimize this effect by choosing the lattice photon frequency ui so 
that the ground and Rydberg states have the same frequency-dependent polarizability a{ijj). The 
results are presented for the case of Rb. 

PACS numbers: 03.67.Lx, 32.10.Dk, 32.80.Rm 



Recently, a number of schemes for quantum computa- 
tion with neutral atoms have been proposed ||^, ^, |[ ^, 
1' 1' 0' 1' ^- those schemes, qubits are realized as in- 
ternal states of neutral atoms trapped in optical lattices 
or magnetic microtraps. The two-qubit quantum gates 
are realized using controlled cold coUisions [0, ||, ||, D, 
controlled dipole-dipole interactions ||, |[ ^, or by con- 
ditional excitations of atoms into the Rydberg states by 
a series of laser pulses |^ . This approach to quantum 
computation has many advantages, such as scalability, 
possible massive parallelism, long decoherence times of 
the internal states of the atoms, flexibility in control- 
ling atomic interactions, and well-developed experimen- 
tal techniques. As extensive experimental studies of the 
feasibility of using neutral atoms for quantum computa- 
tion are under way, more detailed theoretical investiga- 
tion of the various schemes is needed. In this work, we 
investigate optimization of the quantum gate scheme re- 
alized by excitations to Rydberg states We refer to 
this scheme below as the Rydberg gate. The choice of 
this particular scheme results from its potential for fast 
(sub- microsecond) gate operations. 

In the Rydberg gate scheme, the basic qubit is based 
on two ground hyperfine states of neutral atoms confined 
in an optical lattice. A two-qubit phase gate may be re- 
alized by conditionally exciting two atoms to low-lying 
Rydberg states. Different versions of the scheme have 
been proposed In those schemes, either one or both 
atoms may occupy the Rydberg state for much of the 
duration of the gate operation. However, an atom in a 
Rydberg state will, in general, move in different opti- 
cal lattice potential than that experienced by the ground 
state. Therefore, the vibrational state of the atom in 
the lattice may change after the gate operation is com- 
pleted, leading to decoherence due to motional heating. 
The optical potential for a given state depends on its ac 
polarizability, so we can seek to minimize this motional 
heating effect by the choice of a particular Rydberg state 
or of the lattice photon frequency, u. In this paper, we 
describe a method for accomplishing this by matching the 
frequency-dependent polarizabilities, a(w), of the atomic 
ground state and Rydberg state. The results are pre- 



sented for case of Rb; however, the approach used here is 
applicable for logic gates with other alkali-metal atoms. 

In this work, we calculate polarizabilities of the Rb 
atom in its ground state, a^si^), and in various Rydberg 
ns states ans^uj). We demonstrate that there exists a 
lattice frequency lu such that 



a5s(w) ~ a„s{uj). 



(1) 



This finding relies upon the fact |To| that the frequency- 
dependent polarizabilities of alkali Rydberg states are 
close to that of a free electron, afreo(i^) = — e^/moW^, 
where e is the elementary charge in statcoulombs, and me 
is the electron mass. For the optical frequencies treated 
in this paper, ani{i-u) for Rydberg states is comparable 
in magnitude to a5s(0), but it has a negative sign. As a 
function of lu, a^si'-^) increases as co increases from zero, 
but changes sign when uj exceeds the first resonance fre- 
quency. We show that it will always be possible to find a 
frequency for which the Rydberg and ground-state polar- 
izabilities are equal. An optical lattice constructed with 
light of that frequency will have well depths characteristic 
of far-off-resonance traps. 

An alternative approach to polarizability matching of 
ground and Rydberg states is to find Rydberg states with 
resonant frequencies close to those of the ground state. 
In such cases the polarizabilities can be matched at near- 
resonant frequencies, resulting in a tighter trap for a 
given laser intensity. Such possibilities can be identified 
by a straightforward search of tabulated atomic energy 
levels, and they turn out to be rare. The two cases of 
interest in Rb are discussed at the end of this paper. 

We begin with a description of our calculation of the 
polarizabilities. It is convenient for this purpose to use 
the system of atomic units, a.u., in which e, rrie, 47reo and 
the reduced Planck constant h have the numerical value 
1. Polarizability in a.u. has the dimensions of volume, 
and its numerical values presented here are thus mea- 
sured in units of Og, where cq « 0.052918 nm is Bohr 
radius. The atomic units for a can be be converted to SI 
units via a/h [Hz/(V/m)2]=2.48832x 10"^a [a.u.], where 
the conversion coefficient is 47reoag//i and Planck con- 
stant h is factored out. The atomic unit of frequency uj 
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is Eh/h, « 4.1341 X 10^^ Hz, where Eh is Hartree energy. 

The valence contribution to the dynamic polarizabiUty 
for an alkah atom in an ns state can be calculated using 
the formula (in a.u.) jlj] 



1^1 (-g»'pi/, - E,,,){n'p^/2\\D\\ns)'^ 



+ 



„, ^ Ens) 

{^n'p3/2 - Ens){n'p3/2\\D\\ns)^ 



{En 



P'i/2 



E V 



(2) 



where D is the electric dipole operator and Ei is the en- 
ergy of the state i. In this formula, u is assumed to be 
at least several linewidths off resonance with the corre- 
sponding transition. This condition is satisfied for the 
frequencies considered in this work, as the lattice has 
to be significantly detuned from the resonance to mini- 
mize spontaneous emission. The core contribution to the 
polarizability, calculated in the Hartree-Fock (HF) ap- 
proximation, is found to be small (9.3 Og) and is weakly 
dependent on lo in the frequency range considered here. 
The static value for the polarizability of Rb+ calculated 
in the random phase approximation (RPA) ||ll|, |l^ is 
9.1 Oq, close to the value of 9.0 Oq obtained from by Jo- 
hansson from analysis of the observed term values 
of nonpenetrating Rydberg states. The accuracy of the 
core contribution is estimated to be 5 % in Ref . 1 14 . We 
use the RPA value for the core polarizability of Rb+ as 
a baseline, and adjust it to account for valence electron 
and the frequency dependence by using HF calculations. 
The RPA and HF values differ by only 2 %. 

In calculating the valence contribution, the term with 
n = 5 dominates the the sum over states, and the con- 
tributions from higher values of n converge rapidly for 
all values of w considered in this work. The 5s — 5p ma- 
trix elements were measured in Ref. [|l^ and the 5s — np 
matrix elements with n ~ 6, 7, 8 were calculated using 
single-double all-order method in Rcfs. ||lj, |l^. We have 
conducted the calculation of the ground (5s) state a^si^) 
in Rb, including terms with n = 5, 6, 7, and 8 using ma- 
trix elements from Refs. |l5| , [l^ and experimental en- 
ergies from |l^. The contributions to the ground state 
polarizability from states with n' ~ 8 — oo, including 
the continuum, and from terms with n — 2,3,4 are es- 
timated to be very small (0.2 Oq and -0.3 Aq, respec- 
tively). The results for the values of oj near first two 
resonances, 5s — 5pi/2 and 5s — 5p3/2 are listed in Ta- 
ble |. As expected, the value of the frequency-dependent 
polarizability varies very rapidly in the vicinity of the 
5s — 5pi/2 and 5s — 5p3/2 resonances. The second set 
of values of a5s(w) listed in column four illustrates the 
existence of the a^si^) = point between the 5s — 5pi/2 
and 5s — 5p3/2 resonances. The next resonance occurs at 



LU = E; 



5s 



E, 



6Pl/2 



Molof et al. measured a^siO) to be 47.3 ± 0.9 A^, 
equivalent to 319 ±6 Oq. Our calculated value agrees with 
this result to within the experimental uncertainty. The 
static polarizability of the ground state of Rb, 05^(0), 



TABLE I: Dynamic polarizabilities a{Lu) of the ground state 
of Rb in units of ajj; lo is in a.u. (Eh/h). 



LU 


a 


LU 


a 


LU 







318.5(6) 


0.0576500 


-815 (23) 


0.0578 


3931(20) 


0.02000 


360.0(7) 


0.0576645 


-290 (23) 


0.058 


10755(31) 


0.04000 


597.8(8) 


0.0576700 


-97 (22) 


0.059 


-11548(16) 


0.04298 


693.5(9) 


0.0576722 


-21 (22) 


0.060 


-4737(5) 


0.05000 


1211(1) 


0.0576728 


(22) 


0.065 


-1206(1) 


0.05500 


3132(3) 


0.0576734 


21 (22) 


0.070 


-667.1(9) 


0.05700 13854(23) 


0.0576800 


246 (22) 


0.080 


-330.3(8) 


0.05750 


-9311(57) 


0.0576900 


581 (22) 


0.090 


-205.9(8) 


0.05760 


-2869(29) 


0.0577000 


907 (21) 


0.100 


-138(1) 



was previously calculated using the techniques described 
above in Ref. [|ll] . The general behavior of the frequency- 
dependent polarizability of the ground state of Rb has 
been investigated previously in A result for one 

particular value of w = 0.04298 a.u., corresponding to 
A — 1.06/im |2^, has been calculated in the same work, 
using a model potential method. We compare our value 
of the ground state dynamic polarizability for this fre- 
quency with the result of Marinescu et al. [|l^. Our 
resuh a = 693.5(9) ag, hsted in Table |, IS m good 
agreement with the result a = 711.4 Qq from Ref. [ p^ . 
It is just outside the range of uncertainty of the value 
a — 769 ± 61 Aq inferred by Bonin and Kadar-Kallen 
from an atomic deflection experiment ||2l|| . A subsequent 
publication 1 2^ will provide further details of the polariz- 
ability calculations, as well as calculation of the radiative 
widths of np levels to quantify the behavior of the ground 
state polarizability near the resonances. 

Second, we calculate the dynamic polarizabilities of the 
Rydberg states for a specific value of the lattice photon 
frequency lu and then investigate the dependence of po- 
larizability on LU for several Rydberg states. The most 
likely realization of the Rydberg gate scheme involves 
two-photon transitions from the ground state to either 
ns or nd states. For clarity, we calculate the polariz- 
abilities of the ns states in this work. To make an es- 
timate of the ns Rydberg state polarizability ans(j-^) we 
first calculate relevant matrix elements and energies in 
Hartree-Fock approximation. The HF calculations are 
done on a non-linear grid, and the question of the nu- 
merical stability of the calculation of the properties of 
the Rydberg states was investigated by conducting the 
calculation with several different grids to insure the sta- 
bility of the results. We use the resulting matrix elements 
to calculate frequency-dependent polarizabilities of the 
ns states in HF approximation. The results for the ns 
state polarizabilities with lu — 0.0576645 a.u. are listed 
in Table || in the columns labeled q;hf- This frequency 
corresponds to the value of the 05^(0;) = —290 Oq and 
to a detuning approximately 1/3 of the distance between 
the two resonances; aircci^) = —301 Oq at this frequency. 

We observe that there is no significant change of the 
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TABLE II: Dynamic polarizabilities a„s(tj) (in units of ag) 
for Rb, cj=0.0576645 a.u., i.e. A = 790 nm. 



N 



9 
10 
11 
12 
13 



-304 
-292 
-289 
-288 
-287 
-287 



-295" 



14 
15 
16 
17 
18 
19 



-286 
-285 
-284 
-282 
-280 
-277 



"High-accuracy value obtained using experimental energies and 
all-order matrix elements for the dominant terms with n' = 7, 8. 



polarizability values after n > 8. The summation over 
n in Eq. (y) is truncated at n = 23. The contribu- 
tion of the states with n > 23 and continuum is eval- 
uated by carrying out a calculation of the polarizability 
of the 10s state with B-spline basis set The summa- 
tion over the entire basis set yields the result —292 a.u. 
which differs from the value in Table || by only 1 %. To 
evaluate the uncertainty of the HF approximation fur- 
ther we repeat the calculations using HF matrix elements 
and experimental energies from Ref. in Eq. (|). We 
find a substantial, from 10 % to 30 %, difference be- 
tween these two approximations despite only a few per- 
cent differences between Hartree-Fock and experimental 
[Ens — Enp] energies. Such a large discrepancy is ex- 
plained by severe cancellations of the different terms in 
the sum of Eq. (J^ . As an illustration, we list the contri- 
butions to the polarizability of the 15s state from several 
dominant terms together with the values for the corre- 
sponding dipole matrix elements (np||D||15s), energy dif- 
ferences SE — Enp — Ei^s, and denominators of Eq. ^) in 



Table III . We find that the dominant contributions come 
from n — 14: and rt = 15 terms which have different sign 
owing to the sign change in the energy differences. We 
also note that the denominator {SE)'^ — lo^ is completely 
dominated by the term, and, therefore, nearly identi- 
cal for all of the contributions. We also note the cancel- 
lations of the smaller pairs of terms with n = 13, 16 and 
n = 12, 17. 

The cancellations of dominant contributions have also 
been observed in the calculations of the vector tran- 
sition polarizability P and parity-nonconserving ampli- 
tude in Cs and Fr ( see, for example |lj, ^ and ref- 
erences therein) when calculated using similar "direct- 
summation" method. Generally, the values obtained us- 
ing sets of data of the consistent accuracy are more reli- 
able in such cases, so we use HF data below. The accu- 
racy of the values can be increased by estimating corre- 
lation correction contributions to the values of the dipole 
matrix elements. We conduct such calculation for a^si^) 
for the same frequency as data in Table |l| using exper- 
imental energies and high-accuracy all-order matrix ele- 
ments for the dominant terms with n' = 7, 8. The result 
(-295 Oq) is listed in Table || and agrees very well with 



TABLE III: Contributions to the dynamic polarizability of the 
15s state of Rb in a.u., w=0. 0576645 a.u., SE = Enp — Ei^a, 
D = (np||D||15s). The contributions from different terms are 
given in column a(contr). The accumulated values are given 
in column Qf(acc). 



np 


D 


SE 


[SEf -uj'^ 


a(contr) 


a(acc) 


12pi/2 


-5.9 


-0.00215 


-0.00332 


8 


28 


12p3/2 


-8.2 


-0.00214 


-0.00332 


14 


42 


13pi/2 


-15.6 


-0.00111 


-0.00332 


27 


69 


13p3/2 


-21.7 


-0.00110 


-0.00332 


52 


122 


14pi/2 


-111.8 


-0.00034 


-0.00333 


424 


545 


14p3/2 


-161.0 


-0.00033 


-0.00333 


858 


1403 


15pi/2 


144.0 


0.00026 


-0.00333 


-536 


867 


15p3/2 


201.2 


0.00026 


-0.00333 


-1072 


-204 


16pi/2 


16.1 


0.00073 


-0.00332 


-19 


-223 


16p3/2 


23.9 


0.00073 


-0.00332 


-42 


-265 


17pi/2 


6.4 


0.00110 


-0.00332 


-5 


-270 


17p3/2 


9.7 


0.00110 


-0.00332 


-10 


-280 
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FIG. 1: Dynamic polarizability q(lj) for the ns states of Rb 
in atomic units. 



Hartree-Fock value. We note that substituting HF en- 
ergy by experimental values and leaving matrix elements 
unchanged yields substantially different (-356 Qq) value, 
thus confirming our conclusion that using HF values for 
both matrix elements and energies produces more accu- 
rate values then replacing HF energies by experimental 
results. We note that current accuracy of a„s(w) (esti- 
mated at 10 %) is sufficient for the purpose of present 
paper since the polarizability of the ground state varies 
vary rapidly with uj as illustrated in Table | and the vari- 
ation of the value lo within the uncertainty of the Ryd- 
berg state polarizability is small. Consequently, the exact 
point where a^s^Lo) = a„s(tj) can be determined exper- 
imentally by detuning of the lattice frequency near the 
matching point. Further investigations of the Rydberg 
state polarizabilities are beyond the scope of this paper 
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FIG. 2: Dynamic polarizability a{uL>) for the ground state of 
Rb in a.u. The horizontal dashed hne corresponds to polariz- 
abihty of the Rydberg ns level. 



FIG. 3: Levels with transitions energies near 5s — 5p3/2 
resonance. Levels are not to scale. 



and will be considered elsewhere p2[ . 

Finally, we investigate the dependence of the Rydberg 
state polarizabilities on the value of to in the vicinity of 
the 5s — 5pi/2 and 5s — 5p3/2 resonances. As noted above, 
the polarizability of ns state with n > 8 does not change 
significantly with n. The results of the calculation of 
valence dynamic polarizability for the 8s, 15s, and 20s 
states with uj — 0.056 — 0.060 a.u. arc illustrated in Fig. |^. 
We find that the dynamic polarizabilities of the Rydberg 
ns states vary very weakly with uj. 

Summarizing the results above we find that the values 
of the ground state and Rydberg ns state polarizabili- 
ties can be matched near the point between 5s — 5pi/2 
and 5s — 5p3/2 resonances where the ground state po- 
larizability changes sign as demonstrated in Fig. |^. The 
horizontal dashed line corresponds to polarizability of the 
Rydberg ns level. The resonances are shown by the ver- 
tical dashed lines for clarity as we assumed to be a 
few linewidths from the resonances. We predict that 
value of CO where these polarizabilities match equals to 
Auac = 790.14(2) nm for 15s state. This value includes 
10 % uncertainty in the value of the 15s state polarizabil- 
ity and 6 % uncertainty in the value of the ground state 
polarizability. The next point at which the polarizabili- 
ties can be matched is close to the 5s — 6pi/2 resonance 
which requires much larger detuning. The above discus- 
sion is valid for higher Rydberg states, which may be 
chosen for the gate implementation because of the longer 
lifetimes ^J], since we found Rydberg state ac polariz- 
ability to vary weakly with the principal quantum num- 
ber. 

Next, we consider an alternative approach to the is- 
sue of the matching polarizabilities of the ground and 
arbitrary Rydberg state. This discussion is not limited 
to ns Rydberg states. It follows from the expression for 
the dynamic polarizability that if the energy difference 



Ens ~ En'p or End ~ -E'n'p IS accidentally close to the value 
of either 5s — 5pi/2 or 5s — 5p3/2 transition energy then 
the polarizability of the corresponding s, p, ot d state can 
be made large by detuning to the appropriate frequency. 
We have investigated the spectrum of Rb to locate energy 
differences which are close to the resonance 5s — 5pi/2 or 
5s — 5p3/2 transition energies. The closest matches are 
6s — 15p and 4d — lip transition energies which differ 
from the 6E = E^s — E^p^^^ by 20 cm~^ and 60 cm^^, 
respectively. The corresponding energy level scheme is 
illustrated in Fig. |3[ However, in both of these cases the 
terms with such denominators contribute to the polariz- 
abilities of the p states. The required detuning in this 
case is smaller but only lip and I5p levels may be used 
in this method. We note that there are no such Na levels 
for the situation where the lattice light is near Na(3p) 
levels. The advantage of this scheme is high value of the 
ac polarizabilities at the matching point; in the previous 
scheme the value of the ac polarizability at the match- 
ing point is relatively small, leading to the higher laser 
power requirement and subsequent higher scattering rate 
(on the order of 10 s~^). The disadvantage of using np 
vs. ns states in a Rydberg gate scheme is that excita- 
tion of transitions from the ground state requires either 
ultraviolet radiation or a three-photon process. Another 
difficulty of this approach is the necessity of the very fast 
gate operation times to achieve high gate fidelity because 
of the short, 4 fis lifetime of the 15p level. 

We note that by choosing the appropriate Rydberg 
level and a longer wavelength trapping laser (e.d. CO2 
laser) it should often be possible to find a matching 
scheme. In these situations, because of the number of 
nearby resonances for the Rydberg state, one probably 
would have to consider contributions from several inter- 
mediate states. 

In summary, we have found two ways to match 
frequency-dependent polarizabilities of the Rydberg state 
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and the ground state for optical lattices tuned near the 
Rb(5p) states. First, the value of the ground state po- 
larizability at X^ac = 790.032(8) nm crosses zero and 
changes sign. Therefore, it can be matched with polar- 
izabilities of various Rydberg levels near that point as 
demonstrated in Fig. ||. For the case of 15s state, po- 
larizabilities match at X^ac — 790.14(2) nm. Second, the 
value of to can be detuned to nearly match the values 
of the 15p — 6s and lip — Ad transition energies which 
are accidentally close to the 6E = E{5p3/2) — E{5s) en- 
ergy. Then, the dynamic polarizability of the 15p or lip 
state can be made large enough to match the polarizabil- 
ity of the ground state near 6^3/2 — 5s resonance. The 
matching of the frequency-dependent polarizabilities of 
the atom in its ground and Rydberg states results in the 
matching of the optical potential seen by the atom during 
the gate operation and provides the optimal scheme for 
the Rydberg gate operation with respect to the motional 
decoherence. 

Another way of eliminating the differential in the trap- 
ping potential between the ground and Rydberg states is 
to switch the trap off during the time of the gate ac- 
tion, and turn it back on after the action is completed. 
This procedure also induces heating of the atomic center 
of mass motion, the effect of which can be character- 
ized as follows. Suppose that, at time i = 0, an atom 
of mass M is in the ground state of a trap site, where 
for convenience we assume the potential to be that of an 
isotropic oscillator of frequency luq, with ground state en- 
ergy Eq = 3huJo/2. If the trap is turned off suddenly, the 



wavefunction ^'{r, t} for the center of mass coordinate r 
evolves during time t as 



*{r,t} = 



exp(-r2/2(ig(l + iuot)) 

7r3/4d3/2(l + ia;oi)3/2 



(3) 



where do — ^JTi/Mloq. In free expansion of the atomic 
wavepacket from the released trap, its mean kinetic en- 
ergy, < T >— Eq/2 remains constant, while its mean- 
square radius grows as < r^(t) >= 3dl{l+ujQt'^)/2. If the 
trap is suddenly turned back on again at time r, the mean 
energy of the wavepacket in the restored trapping poten- 
tial isE = Eo/2 + MujI < r2(r) > /2 = Eo{l + lj^t^ /2). 
This corresponds to a heating of fceT = {hujQ)LLjQT'^ / 4 per 
cycle of trap release and restoration. Such heating can be 
avoided by use of a matched potential scheme of the type 
proposed here. For a trap frequency of a 1 MHz and a 
1 /iS gate time one obtains fc^T « 0.006(?ia;o)- The gate 
operation time in principle can be optimized by appro- 
priate choice of the Rydberg level, the applied dc field, 
and laser power. 
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